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The thermal one- and two-graviton Green's function are computed using a temporal gauge. In 
order to handle the extra poles which are present in the propagator, we employ an ambiguity-free 
technique in the imaginary-time formalism. For temperatures T high compared with the external 
momentum, we obtain the leading T 4 as well as the subleading T 2 and log(T) contributions to the 
graviton self-energy. The gauge fixing independence of the leading T 4 terms as well as the Ward 
identity relating the self-energy with the one-point function are explicitly verified. We also verify 
the 't Hooft identities for the subleading T 2 terms and show that the logarithmic part has the 
same structure as the residue of the ultraviolet pole of the zero temperature graviton self-energy. 
We explicitly compute the extra terms generated by the prescription poles and verify that they do 
C**) ■ not change the behavior of the leading and sub-leading contributions from the hard thermal loop 

' region. We discuss the modification of the solutions of the dispersion relations in the graviton plasma 

, induced by the subleading T 2 contributions. 
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One of the main motivations for the first attempts to compute the self-energy at finite temperature was the study of 
s ! ■ dispersion relations of a graviton plasma and the related interesting phenomena of anti-damping and wave propagation 
For temperatures T high compared with the external momentum, but well below the Planck scale, the complete 
tensor structure of the leading one- loop contributions, proportional to T 4 , was calculated for the first time in Ref. 0. 
Later some subleading contributions of order T 2 were computed, including the contributions of thermal scalar matter 
and radiation and subsequently all terms proportional to T 4 , T 2 and log(T) where computed taking into account 
' thermal loops of gravitons jj] . 
£N| \ When the internal graviton lines are included, the gauge dependence which arises from the choice of gauge fixing in 
the gravitational action becomes an issue. In Ref. Q the graviton self-energy was computed employing the Feynman- 
de Donder gauge with an arbitrary gauge fixing parameter. While the subleading contributions are gauge dependent, 
the leading T 4 contributions to the self-energy as well as to the one-point function are gauge fixing independent and 
O ' satisfy the Ward identity. This last property is also true for the contributions from matter and radiation, being 
fD , consistent with a gauge invariant effective action for hard thermal loops interacting with gravity. 

• One can go further into the question of gauge dependence by considering a class of non-covariant gauges of the 
£a! \ kind that has been employed in gravity at zero temperature 0, IS IS El IH3 ■ At finite temperature non-covariant 
■ temporal gauges would be even more appropriate, since the Lorentz covariance is already broken by the heat bath 
' but the rotational invariance is preserved. Despite the other well known advantages of the temporal gauge, finite 
temperature calculations have been performed only in Yang-Mills theories both in imaginary and in the real time 
formalisms 0, 0, 0, Q, 0, 0, 0] . This can be partially understood in view of the complexity of the gravitational 
interaction and so explicit calculations in non-covariant gauges have been restricted to the zero temperature case. 
Another reason for the lack of popularity of the temporal gauge in gravity is that, in contrast with the situation in 
Yang-Mills theory, the zero temperature graviton self-energy is not transverse 0,I3- However, this should not be a 
very important concern in the finite temperature case where the transversality property is expected to be violated 
in general. A more important difficulty in the temporal gauge is the problem of spurious singularities arising from 
the n = terms in the Matsubara sums |l8j| . which is even more severe in the case of gravity in view of the higher 
powers of n in the denominator of the temporal gauge graviton propagator. This situation was improved after the 
develo pme nt of an ambiguity-free technique to perform perturbative calculations at finite temperature in the temporal 
gauge |l5l |T(j . Originally this technique was tested using zeta- functions to compute the Matsubara sums and later it 
was applied to the calculation of the gluon self-energy using the standard method of introducing thermal distributions 
by replacing the Matsubara frequency sum with a contour integral in the complex plane of the zero component of the 
internal momentum [T^ j. 

The purpose of the present work is to apply the Lcibbrandt's prescription to the calculation of the thermal one- 
and two-graviton Green's functions in a class of temporal gauges. We will show explicitly how this approach leads to 
a well defined result which can be expressed in terms of forward scattering amplitudes of thermal gravitons |2(i| plus 
contributions from prescription poles. We will also show how the ghost interactions effectively decouple leaving only 
thermal gravitons in the forward scattering amplitudes. We provide the explicit results for the leading and sub-leading 
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hard thermal loop contributions and show that the prescription poles do not change the hard thermal loop behavior. 

In section II we will present the Lagrangian for the graviton field and the corresponding Feynman rules in a class 
of temporal gauges. We will also illustrate the basic approach with the simplest one-loop calculation, namely the one- 
graviton function (tadpole). In section III we describe how the thermal graviton self-energy can be split in two parts. 
The first part arises from the on-shell poles of thermal graviton, and it is expressed in terms of forward scattering 
amplitudes, while the second part is generated by poles in the complex energy plane which are characteristic of the 
temporal gauge prescription. We obtain from the forward scattering amplitudes the leading T A and the sublcading 
T 2 and logT contributions. In section IV we explicitly calculate the contributions from the prescription poles and 
compare the results with the high temperature limit of the forward scattering expression. In section V we employ 
the hard thermal loop results, up to the subleading T 2 contributions, to investigate the modification of the solutions 
of the dispersion relations in a gravitational plasma. We will discuss our results in section VI. 



II. LAGRANGIAN, FEYNMAN RULES AND BASIC DEFINITIONS 

The graviton field, (j)^, can be defined as a small perturbation around the flat space-time metric, rj^ v , as follows 

9/iv (x) — rj^u + K<j>p V (x), k 2 = 32itG. (2.1) 
Here G is Newton's constant and g M „ is the metric tensor. The Einstein Lagrangian density is given by 



C=^yf=ggTR lu , (2.2) 



where is the Ricci tensor given by 



p n pa _ a pa pa p/3 . pa p/3 

r£„ = ^(d^g^ + dvgp.-dpg^) (2.3) 

It is clear from the previous expressions that the Einstein Lagrangian is an infinity series in powers of k (an infinity 
number of terms arises both from the inverse metric g^ v and from the determinant g). Each power n n will come 
out multiplied by a combination of tensor scalar products of n tensor fields <j) and two derivatives dip. Performing a 
systematic expansion in powers of the coupling constant k, it is straightforward to obtain the tree-level Feynman rules 
corresponding to the terms which are quadratic, cubic, etc |3Cj . Before we show the explicit form of these vertices, let 
us recall that the invariance of the Einstein action under general coordinate transformations (gauge transformations) 
imply the existence of Ward identities relating all the vertices down to the quadratic term (see the appendix Ia"|> . 
The identity given by Eq. lAlOfl shows explicitly the usual problem of inverting the free quadratic part of a gauge 
invariant Lagrangian. Following the standard procedure of introducing a gauge fixing condition and ghost fields, we 
add the following two terms to the Einstein Lagrangian |2lJ 

C fix = —tTin^in^) (2.4) 



and 



^-nYf^, (2.5) 

where x M and are the gravitational Faddeev-Popov ghost vector fields. n M is the axial vector and a is a constant 
gauge parameter. Using Eq. (|A3|) . we obtain the following explicitly form for the ghost Lagrangian 

Cghost = X v [n\ d v + t)\ v n ■ d + k (n M (p^dy +<j>„\n-d + ?i M (d\^u))] V X - ( 2 -6) 

Notice that, unlike Yang-Mills theory, ghosts remain coupled to the gravitons even for the choice a = 0. However, 
our explicit calculation will show that the decoupling occurs when the loop integrations are performed. 

We have now all the basic ingredients to perform pcrturbative calculations in thermal gravity. The graviton 
propagator can now be obtained inverting the quadratic term of C + C fi X . Our choice of gauge fixing is such that 
even the bare graviton propagator is already dependent on fourteen independent tensors as shown in table 1 (we will 
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employ this same basis in order to obtain the tensor structure of the thermal self-energy) [3fJ ■ Using this tensor basis 
it is possible to obtain the following compact form for the graviton propagator 
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(2.7) 



where 



Z. 1 
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are convenient linear combinations of the tensors in Table I. As we can see the graviton propagator has the usual 
poles at k 2 = as well as the poles at k ■ u = k = 0. 
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TABLE I: The fourteen independent tensors built from r\p V , k^ and = n M /no and satisfying the symmetry conditions 



pv,pa *vp,pa ^p,v,ap per, pv 



The first and second order terms in k yield the following three and four graviton vertices respectively 
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FIG. 1: Diagrams contributing to the one-graviton function in the one-loop approximation. The curly lines represent gravitons 
and the dashed lines represent ghosts. 



+2k 2a (k 3/3 (r]x^r]ps - VpWSj) + 2k 3p (npxmy ~ ^psmi)) 
+2k 2p (2k 3a ri l 3xr]Sj + hx^VPyVaS - VapVSy)) + 2k 2S k 3p 
x(t7q/3T?a 7 - 2r)pxVa-f)} + symmet. on (a «-> 0), (p <-► A), 
, (S <-> 7)} +permut. of (ki,a,/3), (k 2 ,p, A), (k 3 ,S,j), (2.8) 



V a0,pX,Sy,Ta( k li k 2,k 3 ,k 4 )= — { [k 3 ■ fc 4 ( (VapVpX ~ 2-q ap r)px) {VS^Vtc ~ VSrVya) + 8(V<x8V0p 

+VapVi3S - rj a ^r]sp){r}x n r} T(X - rj^jTjrx) + SVprVasiv^VcrX - Vp^hx)) 
+4:k 3a (2k 4p r]0x ~ kipr) p \)(r) Sl r) Ta - n &T r\ ia ) + 16(k 3p k 4 

aV(38 k 3a k4pT)Sp) 
xiVlvVXr - VlWra) + 8{k 3a k 4s + k 3S k 4a ) (^^7 ~ 2l] 1 xV0p)Vra 

+16k 3a k 4S (r]p T (2ri l3a 'q 7 x ~ Vp-yVvX) + Vi<y( 2 V P rVi3X ~ VpMPt)) 
-16fc 3(5 fc4 Q 77 pr ?7 ( 3 7 ?7 (T A + 2(k 3r k 4& r\ lcr - k 35 ki 1 r} T < 7 )(rio l p'n p x - 2i] ap n f3 x) 
+8(k 3T k 4S r/ 7 x " k 3S k 4:7 rixT){2ri l 3 a r] ap - r) pa r) a p)] 
+symmet. on (a <-> /3), (p «-> A), ((5 <-> 7), (r <-> <r)} 

+permut. of (hi, a, /3), (fo, p, A), (fc 3 , (5, 7), (k 4 , to) . (2.9) 



We have verified that these vertices are in agreement with the Ward identities described in the Appendix El 
Finally, the quadratic and the interacting term in the ghost Lagrangian 1)2.6(1 yields the ghost propagator 



vf: s \k)= l 



kxn p -r] X p, 



2(n.fc) 2 M n.k 



(2.10) 



A// 

and the graviton-ghost-ghost vertex 

V»K, a P u{ki,k 2 , k 3 ) = in{r\ pll n VK n.k 2 + n pli n p k 2L , + n VK n l _ l ki p ) + p <-> k, (2.11) 

respectively. 



A. The one-point function 



In order to introduce our notation and the basic method of calculation we will rederive here the result for the 
thermal one-point function. The one-point function is interesting by itself, since it is directly related to the energy 
momentum tensor derived from the effective action Q ■ It also provides the simplest non-trivial example of a one-loop 
calculation in gravity. Indeed, in contrast with the zero temperature case, the finite temperature one-point function 
is non-zero, being exactly proportional to T A . For that reason it will play an important role in the Ward identities 
obeyed by the hard thermal loop Green's functions. 

The relevant diagrams are shown in the Fig. I. Using the imaginary time formalism PJ] the Eqs. (|2~TU|) and (|2~TT|) 
give the following contribution for the ghost loop diagram shown in figure 1(a) 

r^ ost = kT Y, j f^o _! Vv\ q = 2mnT; n = 0, ±1, ±2, . . . . (2.12) 
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Throughout this work the Matsubara sums like that one in Eq. (|2.12() will be computed using the standard and 
elegant relation [l8j 



-ioo-i-t 



f* [/(«.) + /(-*)] \ coth(^) 



— 200+6 



dgp 
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[/(go) + /(-go)] ( \ + 



CT - 1 



(2-13) 



In general, the vacuum part of the amplitudes (terms which arise from the factor 1/2 inside the round bracket of Eq. 
(I2.13|) ) may be divergent in the limit D — > 4 and so the arbitrary dimension D provides a regulator for the vacuum 
piece of the thermal Green's functions as usual j2^| . 

The tadpole diagram provides the simplest example of an effective decoupling of the ghost graviton interaction 
in the temporal gauge (this is not a trivial property at non-zero temperature). Indeed, substituting (|2.13|) into Eq. 
(|2.12|) we can see that the vacuum piece vanishes as a consequence of the identity 



,D-l 



<fkT = o. 



(2.14) 



The thermal piece also vanishes since we can close the contour in the right hand side of the qo plane without enfolding 
any poles. 

The contribution from the graviton loop in figure 1(b) is a little bit more involved. After some straightforward 
tensor algebra we obtain from Eqs. (|2.7|l and (|2.8|) the following result 



d 



D-l 



q D 



(27T) 



D-l 



2( jD _ 3 )£-|--(D-5)^ 



(2.15) 



It is interesting to notice that the gauge parameter a from the graviton propagator has already canceled out at the 
integrand level. 

Let us now compute Eq. 1)2. 15J) with the help of formula (|2.13(l . As in the case of the ghost loop diagram, the 
contribution proportional to rj^ vanishes. The dimensional regularized vacuum piece will also vanish and we are left 
with only the following expression 



Tui^ — 2 k 



2ixi i 



d J g q^q v 
(2tt) 3 q 2 ' 



Closing the contour in the right hand side plane the pole at qo = |g| gives the following contribution 



o (27r) 3 e -r 



1 



df% q v 



(2.16) 



(2.17) 



where we have introduced g M = qn/\<f\ and J dO is the integration over all directions of q. Finally, using the formula 



dx = T(i>)((v) T v 



we obtain 



7T 2 T 4 



30 



47T 



90 = 191 



7T 2 T 4 



(2.18) 



(2.19) 



where we have employed the quantity u = (1, 0, 0, 0), which coincides with the vector representing the local rest frame 
of the plasma and was introduced in the table I. 



III. THERMAL FORWARD SCATTERING CONTRIBUTIONS TO THE GRAVITON SELF-ENERGY 



The diagrams which contribute to the graviton self-energy are shown in figure 2. The relevant Feynman rules for 
the propagators and vertices are all given in the previous section. Let us first consider the ghost loop diagram shown 
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FIG. 2: Diagrams contributing to the graviton self-energy. The curly lines represent gravitons and the dashed lines represent 
ghosts. The external momentum k is inward. 



in figure 2(a). As we can see from the structure of the ghost propagator in Eq. H2.1U|) the integrand will involve a 
combination of fractions of the following type 



1 



(q ■ u) m [{k + q) ■ u] 7 



0,1,2. 



(3.1) 



Before trying to perform the loop momentum integrations explicitly it is convenient to simplify the integrand using 
well known algebraic identities and change of variables which may reduce the number of terms considerably. Indeed, 
we have found that using a partial fraction decomposition of the quantities shown in Eq. (|3.1|l and a shift q — > q — k 
in the resulting partial fractions containing powers of [(q + k) ■ u]~ , leads to the simplest possible result given by 
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(3-2) 



This procedure (partial fractions and then shifts) has been employed previously in the case of the Yang-Mills theory 
[l7j . In contrast with the present thermal gravity result given by Eq. I)H.2|I . the axial gauge Yang-Mills ghost loop 
vanishes at the integrand level. Notice that the partial fraction decomposition is justified since the integrands are 
regularized accordingly. 

Let us now consider the diagrams shown in figures 2(b) and 2(c). An important difference between these diagrams 
and the ghost loop is that while the ghost loop contains only the poles at go = 0, the structure of the graviton 
propagator in Eq. (|2.7|) is such that there are also the usual simple poles in the right hand side plane located at 
go = |<?| and g = W + k\ — k for the diagram in figure 2(b) and at go = |g| for the diagram in figure 2(c) (notice 
that fco is an imaginary quantity at this stage of the calculation) . In order to use the contour method of integration 
described in section Til Al we will employ the following prescription for the poles at go = [l5l | 
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9S 



<76 



lim 

M-0 (gg - M 2)r 



(3.3) 



With this prescription the temporal gauge poles are moved away from the imaginary axis and we are allowed to 
employ the formula l|2.13|l . The go integral can then be performed closing the contour of integration in the right hand 
side of the go plane, as we did in the previous section in the case of the one point function. The contributions from 
the prescription poles located at go = [i will be analyzed in the next section. 

We now follow the steps explained in the Appendix A of Ref. Basically this consists in the use of Eq. 

(I2.13|l taking into account only the contributions from the poles located on the right hand side plane at go = |g| and 
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<?o = 1 9+ k | — fco . Then, in the residues from the poles at | q+ k \ — kg we perform the shift q — > <f — fc and use the property 
coth(a; + fc ) = coth(x). This yields the following expression in terms of thermal forward scattering amplitudes 



n 



fj,v,a/3\ps 




k, (IV 



90=l«l 



(3.4) 

where the factor 1/2 in front of the curly brackets takes into account the symmetry of the graphs in Figs. 2(b) and 
2(c). It is understood that the external graviton lines with momentum q are contracted with the tensor given by the 
curly bracket of Eq. (|2.7jl . 

We remark that the gauge parameter dependence of Eq. (|3.4|) involves only linear terms in a. This can be 
understood since the quadratic powers of a which could in principle arise from the propagator in Eq. <|2.7[1 do not 
have the on-shell poles. Another interesting property of Eq. (|3.4|1 is that it does not involve thermal ghosts. 

The forward scattering expression in Eq. I|3.4|) is very convenient when considering the hard thermal loop contribu- 
tions which arise from the region where the internal momentum q is of the order of the temperature T, which is large 
compared to the external momentum k. In this regime we can expand the denominators in Eq. i|3.4[) as follows 
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(3.5) 



The leading hard thermal loop contribution is obtained by considering all the integrands which are of degree two in 
the internal momenta q. After some straightforward but very tedious algebra we were able to express the leading 
contribution in the following rather compact form 
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qo=\q\ 



where we have employed the formula (|2.18|l and q have the same meaning as in the Eq. (|2.19|l . 

One can easily verify that this leading T 4 contribution is related to the one-graviton function in Eq. (|2.19|l by the 
Ward identity in Eq. I|A5(I (this result is also in agreement with the calculations performed in the Feynman-de Donder 
gauge @jE1). Since we expect that the leading T 4 contributions are generated by a gauge independent effective action, 
the contributions from the prescription poles in Eq. I|3.3|) should not modify the leading T A behavior. This will be 
confirmed by our explicit calculation in the next section. 

Let us now consider the subleading contributions which are generated when we expand the integrand of Eq. i|3.4|) 
up to terms of degree zero in q. By power counting these will be of order T 2 . In order to obtain the full tensor 
structure generated by the expression l|3.4|l it is convenient to use the following tensor decomposition 
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where the tensors T^'"^ are given in Table I. The coefficients C\ are obtained solving the system of 14 equations 
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where the quantities IT are the following projections of the graviton self-energy 
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Each one of these projections can be expanded using Eq. (|3.5|l . The integrals over the modulus of q can be easily 
performed using Eq. I|2.18|l (they yield the T 2 factor) and the angular integrals are all straightforward. Inserting the 
results for IT into Eq. I|3.8|l and solving for C/, we obtain 
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(3.11) 



There are some properties of the T 2 contributions which are worth stressing. First, the T 2 contributions show their 
gauge dependence explicitly through the gauge parameter a. Each of these gauge parameter dependent terms have 
two powers of momentum relative to the corresponding a-independent ones (the correct mass dimension is provided 
by Uq in the denominators) . Secondly, the simple Ward identity satisfied by the leading T 4 contributions is no longer 
true for the sub-leading contributions. In the Appendix iBl we derive the more general 't Hooft identities and we verify 
that the following identity is satisfied 



A .(0) x\^vaf3 



X 



(3.12) 



where XP/X IS represented by the diagram shown in Fig. 3. The T 2 contribution to xP,\ * s computed in detail in the 
Appendix B. 

We have proceed even further with the hard thermal loop expansion of Eq. (|3.4J) and computed the contributions 
from the integrands of degree minus 2 in q. After integration these yield the logT terms. We have verified that the 
logT contributions of all the projections 11;. (see Eq. (13.9(1 ) arc simply related to the corresponding projections of 
the ultraviolet divergent part of zero temperature graviton self-energy. The zero temperature results were computed 
using the gauge choice a = in Setting a = in our general result we have verified that 
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log 



FS 



log(T)IT 



/3: 



(3.13) 



where II 



is the residue of the ultraviolet divergent zero temperature contribution computed in D = 4 — 2e 
dimensions. The verification of this property in the case of gravity formulated in the temporal gauge complements 
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FIG. 3: The source-ghost diagram. The full/wave line on the left represents the external source. 

similar results obtained in the Feynman-de Donder gauge Q as well as in the case of the Yang-Mills theory [Til . l24l | . 
Since our calculation has been performed for arbitrary values of a, we present complete results in the Appendix C. 

IV. THE CONTRIBUTIONS FROM PRESCRIPTION POLES 



Let us now consider the terms which arise from the poles located at go = M, where fi is the quantity introduced in 
Eq. 1)3. 3|) . It is convenient to express these contributions directly in terms of the projections defined by Eq. (|3.9|) . 
Each one of the fourteen projections can be expressed as follows 
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prcsc 



r=l qo 



q 2 {q + kf 



+ g r l {k ,k- q,q 2 ,k 2 ) 



(go 



(9o 2 -M 2 )'' 



i = 1,2,. 



14. 



(4.1) 



where fT and g\ are polynomials in their arguments, and the denominators <jj have been replaced according to the 
prescription l|3.3|) . This is the most general form of the integrands from the diagrams with gluon or ghost loops. 
Notice that, in particular, the ghost loop expression given by Eq. (|3.2|) yields, after projection, contributions of the 
kind given by the g\ terms above, which contains no on-shcll poles. 

The parameter fi regulates the originally ill defined sums and also makes possible the use of the formula l|2.13[) . 
since now the integrand is regular along the imaginary qo axis. Hence (|4.1I) can be rewritten as 
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2m 
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d D 1 <f— coth. ( ) x 

q 2 \2T)\(ql-^y 
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q 2 (q + k) 2 



(4.2) 



We have employed Eq. I|2.14[) and so the dimensionally regularized integration of the g[ terms has vanished. This 
important property shows how the ghosts arc effectively decoupled at finite temperature. 

Performing the go integration by closing the integration contour at right-hand side plane, we obtain from (|4.2(l 
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In order to obtain the limit /i — > 0, we need the following expansions of the coth and its derivatives 

- +OW , ^ coth (A) 



coth (A) 



d 2 4T 
coth I — i = — 



2TJ 



2T 1 1 
~ + 6T 



<D(n 2 ) 



(4.4) 



coth ( d= ) = - — — + 0{fi 2 ). 



dft 



2TJ 



60 T 3 



An important property of the contribution of the prescription poles is that all the temperature dependence arises only 
from the expansions of the hyperbolic cotangent shown above as odd powers of T . Therefore, the results obtained 
in the previous section for the leading T 4 and sub-leading T 2 and log(T) are not modified by the temporal gauge 
prescription. 

It remains to be verified that the limit fi — > is well defined. Our explicit calculations show that the results for all 
projections do not involve inverse powers of /i. Using the symmetry of the angular integrals all the inverse powers of 
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jj, cancel and we obtain finite results given by 
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(4.5) 



All these integrals are regular and can be done. In the Appendix D we show in an explicit example a closed form 
result. For i — 8, 10, 11, 12 we obtain n^ 1Cf>c = 0, which is in agreement with the 't Hooft identity given by Eq. jB5j). 
We also show that n^ ICSC = for i ~ 7,9, 13, 14. Though the non-vanishing integrals (i = 1, 2, 3, 4, 5, 6) introduce an 
extra temperature dependence, it is clear that they do not change the behavior of the hard thermal loop expressions 
obtained in the previous section. 

We remark that these non-vanishing integrals include both the thermal and the zero temperature contributions 
(notice that the integrands contains a coth instead of the purely thermal part involving the Bose-Einstein distribution). 
Had we computed the thermal part separately we would be left with contributions which are divergent when (i — ► 
as well as the inverse powers of T. Since the dimensional regularization is employed only for the space part of the 
integrals, the vacuum part also contains inverse powers of \x (only the fully dimcnsionally regularized zero temperature 
calculation is well defined in the limit fi — > [l5j]). It is remarkable that the inverse powers of [i in the thermal part 
are exactly canceled by the corresponding ones in the vacuum part and we are left only with the inverse powers of 
the temperature. This property indicates that some of the ill-defined inverse powers of \i have been replaced by a 
thermal regulated expression. In order to understand why these prescription dependent parts are not well defined 
when T — > 0, one should notice that fj, is a dimensionfull parameter which was made "small" in the sense that /.t << T. 
Therefore, the prescription-dependent results cannot be extended to the region T — > 0. 



V. DISPERSION RELATIONS IN A GRAVITON PLASMA 



The sub-leading hard thermal loop contributions proportional to T 2 will produce modifications in the solution 
of the dispersion relations describing the wave propagation in a graviton plasma. The dispersion relations were 
carefully investigated in the case of the leading T 4 contributions Q ■ The inclusion of sub- leading contributions has 
been considered in the Feynman-de Donder gauge Although the sub-leading modification of the solutions of 
the dispersion relations are suppressed by a factor GT 2 <C 1, in relation to the order one part arising from the 
T 4 contributions, one may be interested to know how the gauge dependence of the graviton self-energy will affect 
these solutions (the one-graviton function, which also contributes to the dispersion relations, has no sub-leading 
gauge dependent contributions at the one- loop order considered here). In Yang-Mills theories, the problem of gauge- 
(in)dependence is well understood since a theorem was proved by Kobes, Kunstatter and Rebhan (KKR) [25|]. In 
an one-loop calculation, the gauge dependences of the location of the poles of the gluon propagator are explained in 
terms of the KKR identities. A well known example of this problem is the gauge dependence of the plasmon damping 
constant (see [2(| for a recent review) and its solution by the Braaten and Pisarski resummation scheme (2?J. As far 
as we know, a complete analysis of this problem, in the case of gravity, is still missing. Therefore, we believe that 
it is important to investigate how gauge dependent the graviton propagator is and whether it is possible to extract 
gauge independent information. In this regard, it is remarkable that the one-loop calculations of the QCD damping 
constant, in the axial gauge, though incomplete, satisfy some of the necessary conditions required by any physical 
quantity, being both gauge independent and positive |12j . 

With this motivation, let us apply our axial gauge results in the dispersion relations associated with the transverse 
traceless components of the Jacobi equation for small disturbances in the graviton plasma 0. Proceeding as in 
reference Q, the results given in Eq. I|3.10|l (as well as the corresponding leading T 4 contributions) yields the 
following dispersion relations for the three transverse traceless modes 
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(5.1) 



where L(k) is given by Eq. I|3.11|) . 

Let us now solve these relations in the region of real values of k and k, which is relevant for the propagation of 
waves, and then compare with the corresponding solutions previously obtained in the Feynman-de Donder gauge. It is 
convenient to introduce the dimcnsionlcss quantities k 2 = \k\ 2 /(WirGp), ui 2 = uj 2 / (WirGp) and fig = n^/ (WirGp). We 
will also choose n\ = ui 2 so that the scale of the gauge fixing is compatible with the momentum scale. For intermediate 
values of k and u, the dispersion relations have to be solved numerically and the results are qualitatively similar to 
the ones shown in Fig. 3 of reference 3 . In order to discuss the specific issue of gauge dependence in terms of well 
defined analytic expressions, we will consider the asymptotic regions of very small and very large values of A;, a). In 
the limit fc — > the solution of the dispersion relations l|5.1|) gives the following result for the plasma frequency (this 
is the minimum frequency above which propagating waves are supported by the plasma) 
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where we have neglected higher powers of GT 2 . From Eqs. (4.15) of reference [4j the same limit k — > yields 
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(5.3) 



were ^ is the gauge parameter in the class of covariant gauges and £ = 1 defines the Feynman-de Donder gauge 
employed in reference Q- An important property of these results is that, in both classes of gauges, there is the 
same strong dependence on the gauge parameter. In order to understand this behavior, let us reintroduce the 
dimensionfull parameter WnGp = (8/15) ir 3 GT 4 . Then, in both classes of gauges, one can see that the gauge 
dependent subleading correction is of order (GT 4 )(GT 2 ) which is of the same order as the two-loop corrections, not 
included in this calculation. Therefore, the subleading contributions to the plasma frequency constitute only a partial 
result at the one-loop order. 



In the limit of high frequencies, 



k ^> 1 the asymptotic behavior of the solutions is described by the thermal 



masses mi 



fc 2 (I = A, B, C). Using Eqs. I|5.1() it is straightforward to obtain the following results 
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and 
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When the same derivation is performed using Eqs. (4.15) of reference j4|, for arbitrary values of the gauge parameter 
£, one obtains the following results 
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All these explicit examples clearly show the main differences between these two distinct classes of gauges. It is 
remarkable that the axial gauge subleading contributions contain extra powers of k which makes then larger than the 
corresponding corrections in the covariant gauges, of order (GT 4 )(GT 2 ). Notice however, that the hard thermal loop 
condition, k 2 «T 2 , implies that k 2 GT 2 < 1 so that the subleading contributions will not exceed the leading ones. 
As far as the gauge dependences are concerned, we remark that there are no gauge parameter dependences in the 
axial gauge results (for the choice n 2 , = ui 2 ~ \k\ 2 ). While this property is consistent with the necessary requirement 
that any physical quantity should satisfy, the same is not true when the masses fh\ (I = A, B, C) are computed in the 
covariant gauges. 



In this article we have expli citly computed the thermal one- and two-graviton functions in the temporal gauge. 
We have applied Leibbrandt's |l5j prescription to deal with the temporal gauge poles at finite temperature. This 
calculation provides a rather non-trivial explicit verification of the gauge invariance properties of the hard thermal 
loop contributions. Indeed, the leading T 4 behavior is in agreement with previous calculations in covariant gauges. 
The subleading contributions of order T 2 have a gauge dependence in agreement with the 't Hooft identities. We 
have also compared our logT contributions with the residue of the ultraviolet pole of the dimensionally regularized 
zero temperature graviton self energy, given in Ref. Q, and found that they arc the same (this property has also 
been verified in the Feynman-de Donder gauge Q). Our results include also the full gauge parameter dependence, as 
shown in the Appendix C. 

Our explicit calculation indicates that the temporal gauge may be consistently employed even in the highly non- 
trivial case of thermal gravity. The form of the prescription poles in Eq. (|3.3() do not change the hard thermal loop 
behavior of our main result given by Eq. I|3.4|) . An important property of this forward scattering amplitude is that, 
as opposite to the covariant gauges, it docs not involve thermal ghosts and the gauge parameter dependence is linear 
in a. 

In the analysis of the dispersion relations we have included the hard thermal loop subleading contributions pro- 
portional to T 2 and compared the structure of the gauge dependence with similar calculations which were performed 
earlier in the Feynman gauge. As expected from general formal arguments there are gauge dependent contributions 
which arises from the subleading T 2 terms in the graviton self-energy. By power counting, some of the gauge de- 
pendent terms are of the same order as the two-loop contributions. However, the subleading terms, when computed 
in the axial gauge, are such that their contributions to the asymptotic masses arc enhanced by extra powers of k 
and have a weaker gauge dependence. This behavior is analogous to what happens in QCD, where the plasmon 
damping constant (which also is subleading in the temperature) has a weaker gauge dependence when computed in 
non-covariant gauges. 
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In this appendix we derive the identities which must be satisfied by the vertex functions generated from an action 
which is invariant under coordinate transformations. These identities provide an important consistency check of the 
gravitational Feynman rules as well as the leading high temperature thermal Green's functions. 

The invariance of an action S can be expressed as follows 



VI. DISCUSSION 
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APPENDIX A: WARD IDENTITIES 




(Al) 
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Let us choose the following coordinate transformation with an infinitesimal parameter Se^x) 

x" 1 = x» + fe" (x) . 

Performing the transformation in the metric 

dx a dx^ 

9^u( x> ) = q^-q^jij 9 a p = 9^{x) - g a u(x)d^5e a - g afJ ,(x)d^5e a 

= g^ u (x)+6e x d x g^(x), (A2) 

we obtain 

g'iiA x ) - g^( x ) = Sg^u = k54>^ = -g^a duds' 7 - g va d fi 5e cr - 8e x (d x g f _ LV ) 

= -d v Se„ - d„ Se v - k (d u 8e a ) + <j> va (9 M fe CT ) + 5e x (d x ^ v )] , (A3) 

where we have used Eq. (j2.1(l . Inserting l|A3(l into i|Al|) and using integration by parts, we obtain 

/ d 4 z 5e x (r, vX ^ + r)fi\ d v ) = -« / d 4 x Se x (d„ <f>^ + d» <j> vX + (d x ^)) 4l~ (A4) 
Taking functional derivatives of Eq. \AA\ one obtains the following Ward identities in momentum space 

~X% X (fa)V^ v (fa,fa) = -xU af sx(fa,fa)V 1 ^(fa + k 2 = 0) (A5) 

r\j 

1 



~xl v \(fa)V 3 ^8sJki,k 2 , fa) = -x 1 ^ a 0\(ki,fa)V 2 ^ (k 1 + k 2 ,fa) 

2 fiu 



where 



-xlvSyx(ki,fa)V 2 Zp(h + fa,fa) (A6) 
-xl.u\{kx)V A Zg 5 (fa,fa,fa,fa) = -xlv a i3\{ki,fa)V 3f *g, ' (fa + fa, fa, fa) 

-X^5 7 A( fc 1^3)^ 3 'a/3ra( fc l + ^3,^2, fa) 

-xlvra\(fa^fa)y :i Zl3S 1 ( k i + fa, fa, fa), (A7) 

X%\(fa) = fa^vX+favVnX (A8) 
Xl„a(3\(fa,fa) = favVa^UP + faiiVaWuP + faxVanVPv, (A9) 

and the vertices V n are the momentum space expressions for the n-th functional derivatives computed at M „ = 
(momentum conservation is understood in all identities). 

In the case of a tree-level action there is no one-point "vertex" V 1 ^ and so the quadratic term satisfies the 
transvcrsality condition 

xlux(fa)V^ u (fa,fa) = 0. (A10) 

This may not be the case for an effective gauge invariant Lagrangian. Indeed, it is well know that the one-graviton 
function is non-zero at finite temperature. 

APPENDIX B: GRAVITATIONAL T HOOFT IDENTITIES 

The imaginary time formalism at finite temperature follows closely the corresponding formalism at T = 0. Conse- 
quently, the 't Hooft identities at finite T would be identical to the ones at T = 0, if there were no 1-particle tadpole 
contributions (such terms vanish at T = in the dimensional rcgularization scheme). However, since the tadpole is 
exactly proportional to T 4 , it will not affect the identities involving the sub-leading contributions. To derive these, 
we start from the action 

/ = / d 4 zd 4 y^(x)^ b Q/3 (x - y)4> a p(y) + [ d 4 xd 4 yJ^ (x)X flu x (x - y) V x (y) + ■■■. (Bl) 
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Here denotes the tree order quadratic term plus the sub-leading contributions to the graviton 2-point function 

and X^x represents the tensor generated by a gauge transformation of the graviton field which is given to lowest 
order, in the momentum space, by Eq. (|A8J) . J^ v is an external source, ri x represents the ghost field and • • • stand 
for terms which are not relevant for our purpose. The 't Hooft identity involving the graviton self-energy function is 
a consequence of the BRST invariance of the action /: 

d 4 I t / 7 , -— ^— = (B2) 

In general, Eq. I|B2|) implies the 't Hooft identity 

X^xS^ = 0, (B3) 

which can be written to second order as 

y(°) Y\i^v a[i _ t^(1) T.r2^vaf3 fDA\ 

where V 2lival3 satisfies the identity (|A10j> . Using Eq. (|A10jl we see that l|B4jl leads immediately to the 't Hooft 
identity 

x^ ) A (fc)n s T /3 (fc^)xL%(fc) = o. (B5) 

It is straightforward to show that the identity (|B5|I implies that n su bs = n su bio = n su bn = n su bi2 = and so these 
projections have a temperature behavior which is at most proportional to T 4 . 

In order to verify (|3.12(l we need to calculate the tensor xiuX which appears in l|B4(l . In this way we need the 
source-graviton-ghost vertex which can be obtained from the Lagrangian |28j| 

C s = kJ^D^ (B6) 

Using the transformation Sg^ = nD IJiU \e x we obtain 

/ 



kj 



The diagram in Fig. 3 can be now calculated using vertex l|B7|l and Feynman rules Ij2.7|) . (|2.1()|l and l|2.11|) . 

Expanding XnvX m ^ ne base shown in Table II and using the forward scattering method as we did for the one- and 



k^ky k\ 

k^UyUx + kyU^ux 
u^kykx ~\~ Uyk^kx 
Vv-vkx 
■q^ux 
k^rjvx + ky-q^x 
u^x + u v riv,x 



TABLE II: 10 independent tensors base 



two-graviton functions, we obtain the following leading T 2 contribution for X^lx 



= jg k o [ 7 l\v u n + Vx^u u - 2u f _ l UyU X ] (B8) 



15 



Contracting JB8J with V 2tival3 yields 
,0) -> « 



^JA y2MI/a/3 = ■Y^-^o{2fc-M(fc Q u, 3 + fc / 3U Q ) + A:-u(277 Q(3 A: A -?/ A/ 3fc Q - r/ Xa k p ) 



(B9) 



+ fc 2 (?7qAW,3 + J7/3aw q - 2u a upu\) - k x [k a up + kpu a ) - 2k ■ u 2 rj a pux) 



Using the result for II^ , which can be obtained inserting Eq. (|3.1U|) into (|3.7[) . we have verified that the contraction 



with Xnlx yields Eq. i|B9|) with opposite sign in complete agreement with Eq. (|3.12|> . 



APPENDIX C: THE logT CONTRIBUTIONS FOR ARBITRARY VALUES OF a 

In this appendix we complement the result presented in Eq. (|3.13|) and include the contributions proportional to 
the gauge parameter a. We have obtained the following results for the 14 projections (see Eq. Ij3.9|l ) 
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where we are using the quantity y = -£ in order to compare with the zero temperature results of reference (8j- 

In terms of these projections, the coefficients of the transverse traceless components of the graviton self-energy, as 
defined for instance in reference 0, can be written as 
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These expressions show that the dispersion relations associated with the transverse traceless modes will, in general, 
be gauge dependent at this order of perturbation theory. 



APPENDIX D 



As an example of the calculation shown in subsection 4.2 we will calculate a contribution of the prescription poles 
for the projection IIg rcsc . This contribution can be written as 



n prcsc = K T J2 
fc ° 



dP'^q 1 (£>-3)D 
(2-k) d - 1 2 D-2 



(5 - 2D) fc 4 + (2D - 6) k 2 (fc 2 + q 2 ) + fc 4 - 2g 2 fc 2 
9o (q + k) 2 



(Dl) 
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Using the prescription (|3.3(l and Eq. (|2.13|) we obtain 



T-rpresc 
Xi 6 



^°°+ £ dq f dP~ X q 



ko P^oJ_ ioo+€ Am J (2m) D ~ l \2T ) 2 D - 2 



— lim 



<? 4 + (5 - 2D) k% + (2D - 6) fc 2 , (fc 2 + g 2 ) + fc 4 - 2q 2 k 2 



+ q 



-<i 



(q 2 -^) (q + k) 2 

Closing the integration contour at right-hand side plane and expanding in terms of a power series in fj, we obtain 

d°- V D (D - 3) I (2D - 6) fc 2 (q 2 + fc 2 ) + (5 - 2£>) fc^ + (fc 2 - q 2 ) 2 



(D2) 



T-rprcsc 
11 6 



K 2 T 



(2^)°^ 2 D-2 



k 2 - (q + fc) 2 



fc — fc 



K 2 T / d D ~V D(D-3) 



(2ttz) 



L> - 2 



(2L> - 6) fc| (q 2 + 2fc 2 ) + (5 - 2D) k% + q 4 + 4 fc ■ q 



(kl + q 2 ) 2 



(D3) 



where k± = iko and we have performed a shift q — > q — fc. In the limit fl^4wc obtain 

8C= 2 _£ 2 



(D4) 



The contributions to the projections n^ 1080 , Ilg rc8C , n^3 CSC and n^ 80 are obtained in a similar way and we find that 
they vanish. However, for the projections n^ rci!C , il2 lcsc , rig" 3130 , n^ 10 ^ and ng 1C8C we have more involved expressions 
containing inverse powers of T as in Eq. i|4.5|) . 
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